Force-velocity relations for multiple-molecular-motor transport 
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A transition rate model of cargo transport by N molecular motors is proposed. Under the as- 
sumption of steady state, the force-velocity curve of multi-motor system can be derived from the 
force-velocity curve of single motor. Our work shows, in the case of low load, the velocity of multi- 
motor system can decrease or increase with increasing motor number, which is dependent on the 
single motor force-velocity curve. And most commonly, the velocity decreases. This gives a possible 
explanation to some recent experimental observations. 

PACS numbers: 87.16.Nn 



I. INTRODUCTION 

The cargo transport by single cytoplasmic molecu- 
lar motor has been widely studied both experimentally 
0, 0, S, 0, Hi and theoretically [10, 0. Cargoes in vivo, 
however, are typically transported by several motors M 
and sometimes even by different kinds motors [n| lllj - 
By far little has been known about the cooperativity of 
multiple motors during cargo transport, and it is still an 
important and open research subject. Especially, cargo 
transport by multiple processive motors attracts much 
attention, since these motors can transport cargoes over 
long distances without unbinding from the track, which 
is convenient for experimental studies. Actually, experi- 
mental and theoretical studies on such systems have been 
carried out in the last decade. These works investigated 
systems of fluid-like car goes fl2l . Il3l . [l4j , or rigid or elastic 
cargoes [HI, [TH, [13, EM l2fj |. in the existence of exter- 
nal load force. In the cases of rigid or elastic cargo, the 
cargo can induce strong coupling between motors, which 
is the focus of this article. 

Recently, Klumpp et al. proposed a transition rate 
model to study the cooperative cargo transport by pro- 
cessive motors [21] . In their model, motors are supposed 
to share the load force equally and have no other mutual 
interactions. By this treatment, the author concluded 
that the velocity of the cargo increases with the increas- 
ing motor number. Theoretical analysis based on ratchet 
models also give the same results [III, [H, |2(J. Never- 
theless, experiments have shown that the cargo velocity 
is approximately independent of the number of motors 
pulling the cargo ELLZLSl]. Very recently, Shubeita et 
aUs experiment showed that increase of kinesin number 
leads to slightly reduced cargo velocity (23[. This result 
is out of expectation and contradicts some theoretical 
results PbSl . [l9L [20j , but is supported by the simulation 
results |24| . 

Most of the theoretical studies mentioned above fall 
into special cases since they are dependent on the mod- 
eling of single motor stepping. In this article we want 
to generally investigate the dependence of the velocity of 
multi-motor system on the number of motors pulling the 



cargo. For convenience, we suppose the motors can not 
detach from the track. We proposed a steady state tran- 
sition model of transporting cargo by N motors. Our 
calculation shows that the velocity of N motors trans- 
port depends on the single motor's force- velocity(F-V) 
relation, and especially in the case of low load, the veloc- 
ity of multi-motor system may decrease with increasing 
motor number. This result provides a new explanation 
to Ref. [23| , and our work predicts a general behavior of 
multi- motor transport. 



II. MODEL AND RESULTS 

A. Two-Motor case 

We first consider the situation of two motors (A and B 
in Fig. 1) transporting a cargo. The equilibrium distance 
between A and B is denoted as Iq. If motor A takes a 
forward step, the distance between A and B increases to 
produce a traction interaction between them [25l[26j. and 
the cargo is assumed to step forward with a distance d/2, 
where d is the step size of the motor. This assumption is 
quite reasonable, e.g., it has been shown experimentally 
that in the case of two kinesins transporting a microtube, 
the step size of cargo is 4nm which is half of a single 
motor's step size [l7| ■ When motor B takes a forward 
step, resulting a repulsion interaction between the two 
motors, and the cargo also take a forward step of d/2. 
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FIG. 1: (Color online) A cargo is transported by two motors 
A and B. 
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Here we focus mainly on the cases that the single motor 
performs unidirected stepping, since it's a good enough 
description of processive motors like kinesin. 

States of the two-motor system can be specified by 
the "effective distances" S c g — Sab — h, Sab the real 
distance between the two motors. So each state of the 
system can be denoted by i), where i — S c g/d are in- 
tegers. When motor A takes a forward step, the system 
transits from state |t) to state \i + 1). While motor B 
takes a forward step, the system transits from state \i) 
to state \i — 1). Since the time for either motor perform- 
ing the stepping motions is much shorter than the dwell 
time(e.g., the 8-nm stepping motion of kinesin occurs on 
the microsecond timescale while the dwell time is always 
larger than millisecond 0]. The 36- nm stepping motion 
of myosin-V occurs within a few millisecond, far less than 
the dwell timescale of second the two motors can 

never step forward simultaneously. Therefore, the tran- 
sitions between states can be expressed as 



Now we turn to the mean velocity of the two-motor 
system. Denoting by Pi the probability that the system 
is in state Here we are concerning the steady-state 
velocity of the system. The steady-state solution of the 
process described by Eq. (1) can be expressed as 



Pi = p n 



3=0 W i+1 



for (i > 0), and P_; = P,. (3) 



Considering the normalization y^L Pj = 1, Pq satisfy 



Pn = 



n i— 1 + 



(4) 



In the case of linear spring linkage, when either of the 
two motors takes a forward step, the cargo goes forward 
d/2. So the average velocity of the cargo is then given by 



■••;=!|-i>^|o)^|i)^---^=» 

U3_ x uj u 1 u> 2 uj n 

The minus represents the distance between the two mo- 
tors smaller than Iq. The transition rates tuf 1 between 
states depend on the external load force and the cargo- 
mediated force between the two motors which will be 
discussed below. When stall force is reached on either 
motor, the system gets into the extremity states, \n) or 
|-n). 

The cargo-mediated force exerted on either motor is 
quite intuitive, i.e., when i > 0, motor B exerts a re- 
sistance force / on motor A through the cargo, while 
motor A exerts an assistance force — / on motor B; and 
it is contrary when i < 0. The magnitude of / is deter- 
mined by the distance of the two motors and the stiff- 
ness of motor linkage which connects the motor heads 
to the cargo as shown in Fig. 1. For different kinds of 
motor and cargo, the stiffness of linkage are different. 
Here we take kinesin as an example, and the methods 
can be extended to myosin-V directly, but it does not 
apply to dynein because of their complexity and unclear 
stepping behaviors [H, . The linkage of kinesin exhib- 
ited an adequately linear behavior In such cases, 
the internal force between the two motors can be easily 
expressed as / = idk/2, where k is the linkage stiffness, 
and i times d is the effective distance between the two 
motors and each motor shares one half of the distance. 
When an external load F is taken into account, the to- 
tal force borne by each motor should be (F/2 + f) or 
(F/2 — /). Since the force- velocity relation V\(F) for 
single motor transport has been widely studied both ex- 
perimentally and theoretically, one can easily know the 
step rates R\ (F) — V\ (F) / d for either motor of the sys- 
tem. Then we can get the transition rates tof 1 in Eq. 
(1), 



V 2 (F) 
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[V^F/2 + idk/2) + Vi(F/2 - idk/2)} 



(5) 



Ul; 



R 1 {F/2± idk/2). 



(2) 



Defining V 2 (F) = V 2 (2F). It's obvious that V 2 (F) < 
Vi(F) rigorously holds if the single motor F-V curve 
is purely concave, which is followed by two facts: (1) 
V 2 (F) < Vi(F) when the load F is near zero, i.e, two- 
motor transport is slower than single motor transport 
at low load. (2) V 2 (F) > Vi(F) when F is large, i.e., 
the two-motor transport is generally faster. While single 
motor F-V curve is purely convex, then V 2 (F) > V\(F) 
rigorously holds, and two-motor transport is faster than 
single motor transport in the whole range of F. 

Most real single motor F-V curves, however, are usu- 
ally a mixture of concave and convex regions, so one can't 
intuitively know the characteristics of the 2-motor F-V 
curve from the single motor F-V curve, but can still get 
some insight of V 2 {F) when F is near zero. Roughly 
speaking, we have two typical categories of single motor 
F-V curve. 

Category A: the single motor velocity is more sensitive to 
resisting load than to assisting load (i.e., the F-V region 
of assisting load is concave and much flatter than the 
region of mediate resisting load, as illustrated by Fig. 
2 A. ), V 2 (0) < Vi(0) may usually hold. 
Category B: the single motor velocity is more sensitive 
to assisting load than to resisting load (as illustrated by 
Fig. 2B), V 2 (0) > Vi(0) holds. 

Therefore, for real single motor F-V curves, we can obtain 
similar results as for purely concave and convex curves. 

For the nonlinear-spring motor linkage case, the step 
spacing of the cargo varies. Eq. (5) seems not able to be 
used in this case. But noticing that the average velocity 
of the cargo is equal to the average velocity of either 
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FIG. 2: (Color online) The V N (F/N) curves of N-motor 
system derived from two typical single motor F-V curves 
as shown in insets. The linkage stiffness is taken as k = 
0.3pN/nm just for illustration. The value is adopted from the 
experimental result of Ref. The curve in inset of (A) is 
fitted from the experimental date of Ref. l] with [ATP]=1.6 
mM. The inset of (B) is adapted from theoretical results of 
Ref. with [ATP] =5 fiM. 



motor in the long time limit, so the cargo velocity can be 
expressed as 



y 2 (F) = d^p, 



(6) 

where the transition rates are tuf = V\(F/2 ± fi)/d and 
fi is the internal force between motors for state The 
last term of Eq. (6) is similar to Eq. (5). 

Considering of backward steps, we set the forward step 
rate and backward step rate for a single motor are R{ (F) 
and Rl(F) respectively with Fi(F) = [R{(F) -R\ (F)]d, 
and their ratio is e(F) = R\(F) / R[(F) which has been 
studied in Ref. '2\. If V\(F) and e(F) are given, R((F) 
and R\(F) can be known. The transition rates wf' in 




FIG. 3: Transitions between states for N-motor system. 

Eq.(l) are then 

uf = R{{F/2 ± /,;) + R\{F/2 T fi). (7) 

One gets the probabilities Pi of state \i) by Eq.(3) and 
the average velocity of the two-motor system, 

m j fl Ha±MiUffi:M, , 8 ) 

i— — n 

which is the same form of Eq. (5). 

B. Multi-Motor case 

For a general N-motor system, we can regard this sys- 
tem as the combination of a single motor and a (N — 1)- 
motor subsystem. In order to conveniently describe the 
model, we call this very single motor as motor A, and the 
(N — l)-motor subsystem as "motor" B. In the case of 
linear-spring motor linkage, if one of the N motors take 
a forward step, the N-motor system will progress d/N. 
So the step size of motor A is d, while the step size of 
"motor" B is d/(N — 1). Similar to the two- motor sys- 
tem, we can express the states of the N-motor system by 
the effective distance between motor A and "motor" B. 
Each state denoted as where i is the value of the ef- 
fective distance between motor A and "motor" B divided 
by d/(N-l), i.e. i= {N-l)S/d, where S is the effective 
distance between motor A and "motor" B. When motor 
A takes a forward step, the N-motor system will transit 
from state \i) to \i + N — 1), while "motor" B takes a 
forward step, the system will convert from state \i) to 
\i— 1). Transitions between the states are shown in Fig. 
3. The minus states represent that the effective distance 
between A and B smaller than the equilibrium distance. 
If the effective distance between A and B is S, then this 
distance shared by motor A is S(N — 1)/N and shared 
by "motor" B is S/N, so the internal force between A 
and B is / = idk/N. Then the transition rates between 
the states can be given as 

uj+ = Vi(F/N + idk/N)/d, 

ujr = (N - l)V N -i ({N - l)F/N - idk/N) jd 

= (N- l)V N -i (F/N - idk/(N 2 - AT)) /d, (9) 

where the definition V N -\(F/(N - 1)) = V N -i(F) is 
used. The steady state solution of the transition model 
shown in Fig. 3 can be obtained if V\(F) is given and 
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Vn-i(F) is known by the recursion of Eq. (9) and Eq. 
(10-11) given below. Any one of the motors taking a for- 
ward step will make the system step forward d/N. The 
average velocity of the N-motor system is then given by 

V N (F) = V N (F/N) = ^^Pi K + (10) 

i 

According to the network structure of Fig. 3, we can get 
Ei^r =(N- 1)E,^+ the Eq. (10) can also be 
expressed as 

V N (F) ee V N (F/N) = dJ2^ u+ = ]v3T E Pi U T- 

i i 

(ii) 

If the number of motors N is even, the above descrip- 
tion can be greatly simplified. Dividing these motors 
into two groups, each group contain the same number 
of motors, N/2. Then this N-motor system can be re- 
garded as a two-big-motor system. If the F-V curve of 
the (N/2)-motor system has already been known, we can 
easily obtain the F-V curve for the N-motor system by 
the same method of two-motor system. For convenience, 
we calculate the F-V curves for this situation and show 
the results in Fig. 2. 

Fig. 2A displays the F-V curves of multi-motor trans- 
port, as well as the single motor F-V curve of Category A 
which is the most common case in many previous experi- 
mental and theoretical studies. The calculation shows the 
velocity of multi-motor system decreases with increasing 
motor number in the case of low load. This offers a pos- 
sible explanation to a recent experimental observation 
[HI . We also notice that some experiments conclude that 
the increase of motor number does't affect the transport 
velocity of cargoes [u| [13, This may be a conse- 

quence of the fact that some F-V curves of single motor 
are not far from linear, so all Vn(0)(N — 1, 2 • • • ) are al- 
most equal, i.e., multi- motor transport is not significantly 
slower than single motor transport at low load. 

One can also consider the consequence of motor de- 
tachment. Suppose M motors adhering on a cargo. The 
number N of binding motors is no longer constant but 
varies with time between zero and M. Therefore, the 
mean effective velocity of cargo transported by these M 
motors can be expressed as weighted average of V/v (e.g, 
Eq.(6) in (H), 

M 

y e ff( F ) = J2 v * p n(f) (12) 

N=l 

where Pjf(F) is the force-dependent equilibrium bind- 



ing probability of N motors. One can easily show that 
V^j{F) > V e f 2 f (F) holds at low load if M x < M 2 , by 
noticing that Vn(F) decreases with N at low load as 
shown in Fig. 2A. 

III. DISCUSSIONS 

In this paper the F-V curve for cargo transport by 
multiple motors has been discussed. We focused on the 
linear-spring motor linkage case without considering the 
motor's backward steps. The results show that the F-V 
curves of multi-motor system depend on the F-V curve of 
single motor. Novel insights are gained through our cal- 
culation, i.e., at low load, the velocity of the multi-motor 
transport decreases with the increasing motor number if 
the single motor F-V curve belongs to Category A, and 
increases if the single motor F-V curve belongs to Cat- 
egory B. Our linear-spring motor linkage model can be 
extended to the nonlinear-spring motor linkage case and 
also the case of existence of backward steps. Results of 
both the latter models are qualitatively consistent with 
the result of the former model. 

Our results contradict to earlier results which predict 
multi-motor transport is faster than single motor trans- 
port [HI, EiJ H3|- But the very recent experiment sup- 
ports our results, which shows that increasing of mo- 
tor number causes slight decrease of cargo velocity [23j j . 
Ref. [1H attributes the decreasing of cargo velocity to the 
detachment of motors from filament. From our results, 
even without motor detachment, the cargo velocity can 
still decrease with increasing motor number. In fact, the 
Fig. 2 of Ref. [l6[ also shows a slight decrease of cargo 
velocity with the increasing of motor numbers. There 
is another difference between our results and the results 
of Ref. [24j . In Ref. [24[, the simulation results show 
that the multi-motor transport is slower than single mo- 
tor transport in low load case, but if the motor number 
is larger than two, the cargo velocity will increase with 
the increasing of motor number. Therefore, further ex- 
perimental tests are needed, for example, hopefully by 
the method of Ref. p]}. 
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